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Abstract 

Chemistry-climate models predict an acceleration of the upwelling branch of the 
Brewer-Dobson circulation as a consequence of increasing global surface tempera- 
tures, resulting from elevated levels of atmospheric greenhouse gases. The observed 
5 decrease of ozone in the tropical lower stratosphere during the last decades of the 20th 
century is consistent with the anticipated acceleration of upwelling. However, more re- 
cent satellite observations of ozone reveal that this decrease has unexpectedly stopped 
in the first decade of the 21 st century, challenging the implicit assumption of a continu- 
ous acceleration of tropical upwelling. In this study we use three decades of chemistry- 
io transport-model simulations (1980-2013) to investigate this phenomenon and resolve 
this apparent contradiction. Our model reproduces the observed tropical lower strato- 
sphere ozone record, showing a significant decrease in the early period followed by 
a statistically robust trend-change after 2002. We demonstrate that this trend-change 
is correlated with corresponding changes in the vertical transport and conclude that 
is a hiatus in the acceleration of tropical upwelling occurred during the last decade. 


1 Introduction 

The issue of whether the large-scale Brewer-Dobson Circulation (BDC) has strength- 
ened in the recent past, as a result of anthropogenic activity, has been raised (Oman 
et al., 2009; Butchart et al., 201 0; Randel and Jensen, 201 3). Recent chemistry-climate 
20 model (CCM) simulations predict an increase of resolved wave activity and orographic 
gravity wave drag resulting from increasing sea surface temperatures (Garcia and Ran- 
del, 2008; Oman et al., 2009; Waugh et al., 2009; Butchart et al., 2010; Garny et al., 
201 1 ). This strengthens the upwelling branch of the BDC, commonly referred to as the 
tropical upwelling. In comparison, the behaviour of the observations is ambiguous. The 
25 long-term cooling of the tropical lower stratosphere (LS, about 17-21 km; Thompson 
and Solomon, 2005; Young et al., 2012) and the observed weakening of the strato- 
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spheric quasi-biennial oscillation (QBO; Kawatani and Hamilton, 2013) are consistent 
with the predicted increase of upwelling. On the other hand, the mean residence time 
of air parcels in the stratosphere (age of air) inferred from sulfur hexafluoride (SF 6 ) 
measurements is inconsistent with an overall acceleration of the BDC (Engel et al., 
5 2009; Stiller et al., 2012). They indicate no significant changes or even deceleration of 
the vertical transport in the middle stratosphere. To reconcile the observed discrepan- 
cies it has been argued that the individual branches of the BDC are evolving differently, 
i.e. an increase in tropical upwelling does not necessarily imply an acceleration of the 
overall circulation (Bonisch et al., 2011; Diallo et al., 2012; Lin and Fu, 2013). 
io Ozone (0 3 ) is a sensitive proxy for vertical transport in the tropical LS (Randel et al., 
2006; Waugh et al., 2009; Randel and Thompson, 2011; Polvani and Solomon, 2012). 
Its local mixing ratio is considered to result from a stationary state involving production 
by oxygen ( 0 2 ) photo-dissociation and a steady influx of 0 3 -poor tropospheric air from 
below (Avallone and Prather, 1996; Waugh et al., 2009; Meul et al., 2014). Meridional 
is mixing from higher latitudes is a secondary effect that contributes to the seasonality in 
the 0 3 mixing ratios (Ploeger et al., 2012). Several studies have reported a negative 
trend of 0 3 in the tropical LS in the range of -(3-6)% per decade, consistent with 
the CCM predicted increase of tropical upwelling (Randel and Thompson, 201 1 ; Sioris 
et al., 2014; Bourassa et al., 2014). In contrast, more recent 0 3 observations from 
20 various satellite instruments indicate no statistically significant decrease of LS 0 3 since 
the beginning of the 21st century (Kyrola et al., 2013; Eckert et al., 2014; Gebhardt 
et al., 2014). 

Stimulated by the need to explain the unusual linear trends revealed from the vertical 
profile of 0 3 retrieved from SCI AM ACHY 1 we use three decades of 0 3 observations 
25 and simulations to investigate this phenomenon. Section 2 describes the observations, 
model and regression analysis used in this study. The results are discussed in Sect. 3. 

1 First reported at the Quadrennial Ozone Symposium 2012 Toronto, 27-31 August 2012 
and published in Gebhardt et al. (2014). 
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2 Data and analysis 
2.1 Observations 

For a quantitative analysis of tropical upwelling, we use combined 0 3 observations 
from satellite instruments and sondes. The earlier decades (1985-2005) are covered 
by the ERBS/SAGE II instrument (McCormick et al., 1 989), providing 0 3 profiles based 
on solar occultation measurements. Due to its viewing geometry, the vertical resolu- 
tion of the profiles is high (1 km, range 15-50 km), although the horizontal sampling 
is relatively sparse (global coverage in 1 month). Here we use version 7.0 of the data 
(Damadeo et al., 2013), screened for cloud and aerosol contaminated profiles as sug- 
gested by Wang et al. (2002). Two years of data after June 1991 has been omitted due 
to contamination by the eruption of Mt. Pinatubo. For the last decade (2002-2012), 
we use 0 3 observations from ENVISAT/SCIAMACHY (Burrows et al., 1995) based on 
limb geometry (retrieval version 2.9; Sonkaew et al., 2009). The vertical resolution is 
about 3-4 km over an altitude range of 10-75 km; global coverage is achieved every 6 
days. Data from both instruments has been binned into monthly samples on a uniform 
horizontal and vertical grid (15° Ion. x 5° lat. x 1 km). To minimise sampling issues and 
taking into account the differences in horizontal and vertical resolution of the instru- 
ments, any further analysis is based on partial columns of 0 3 between 17-21 km and 
20° N-20°S, similar to the approach of Randel and Thompson (2011). 

The satellite data is augmented by an ensemble of tropical sonde measurements 
from the Southern Hemisphere Additional Ozonesondes network (SHADOZ; 1998- 
2013; Thompson et al., 2003, 2012). We use 10 sites located in the tropics with long 
and continuous records. The selected stations along with their temporal coverage and 
mean value are listed in Table 1 . Typically there are 2-4 observations per month for 
each SHADOZ station, which provide 0 3 profiles in a considerable higher vertical res- 
olution (50-1 00 m) compared to the satellite instruments. As there is a high degree of 
longitudinal symmetry in the stratospheric ozone profiles (Thompson et al., 2003), we 
average the individual records to obtain a representative mean for the tropics. 
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2.2 Model 

To obtain a consistent timeseries of LS 0 3 of the last decades for direct compar- 
ison with observations, we conducted a 33-year simulation with the Bremen three- 
dimensional chemistry-transport-model (B3DCTM; Sinnhuber et al., 2003; Aschmann 
5 et al., 2009; Aschmann and Sinnhuber, 2013). The current version of the model has 
a horizontal resolution of 3.75° Ion. x2.5°lat. and covers the vertical domain from the 
surface up to approximately 55 km using a hybrid o - 9 coordinate system (e.g., Chip- 
perfield, 2006). The vertical resolution in the tropical LS is about 600 m. The model 
is driven by 6-hourly input of European Centre for Medium-range Weather Forecast 
io (ECMWF) Era-Interim (El; Dee et al., 2011) reanalysis data. Vertical transport in the 
purely isentropic domain (above « 1 6 km in the tropics) is prescribed by El all-sky heat- 
ing rates. The B3DCTM incorporates a comprehensive chemistry scheme originally 
based on the chemistry part of the SLIMCAT model (Chipperfield, 1999), covering all 
relevant photochemical reactions for stratospheric 0 3 chemistry. Reaction rates and 
is absorption cross sections are taken from the Jet Propulsion Laboratory recommenda- 
tions (Sander et al., 201 1 ). To avoid initialisation artefacts, the model has been run with 
replicated input data to reach steady state before starting the actual integration from 
January 1 979 to October 201 3. 

2.3 Regression 

20 The multivariate regression analysis used throughout this study is based on Reinsel 
et al. (2002) with Y x as the monthly mean variable to be fitted: 


Y x — fi + Sj + ® i X-f f + ®2 ^2t QBO( + ENSO t 

+ SC t + A/ t (1) 

25 Here, /i is the baseline constant, S t a seasonal component, &) : 2 are the trend coef- 
ficients with X x 2t as linear trend functions and A/ t represents the unexplained noise. 
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Note that in contrast to previous studies, which examined LS 0 3 (Randel and Thomp- 
son, 201 1 ; Sioris et al., 201 4), our regression model assumes two linear components, 
which account for a possible change of trend. Here, ® 1 is the linear trend up to a spec- 
ified inflexion date T 0 . After T 0 , the new linear trend ® comprises the sum of the earlier 
5 trend and the trend-change component ® 2 . The additional regression terms are 
QBO t for QBO, ENSO, for the El Nino Southern Oscillation (ENSO) and SC, for solar 
cycle. The QBO proxy consists of the QBO.U30 and QBO.U50 (zonal wind 30/50 hPa) 
from the NOAA Climate Prediction Center 2 , the ENSO proxy is represented by the 
Multivariate ENSO Index (MEI) from the NOAA Earth System Research Laboratory 3 
io (Wolter and Timlin, 2011) lagged by two months and the solar cycle by the Bremen 
composite Mg II index 4 (Snow et al., 2014). 

Assuming first order autocorrelation noise (AR(1) model), as commonly used in the 
regression of 0 3 timeseries (e.g., Reinsel et al., 2002; Jones et al., 2009; Sioris et al., 
2014), the corresponding errors for the trend components calculate as 


°N / I + 

^ ”,, 3 / 2 ^ 1-0 

On / l + 0 ( n \ 3/2 

‘ 02 * 2 V 1 - 0 UoflJ 

o N / l +0 ln 0 +4n 1 
^ ~ n 3/2 y 1 - 0 V 4/7 


2 www.cpc. ncep.noaa.gov/data/indices/ 

3 www.esrl. noaa.gov/psd/enso/mei/ 

4 www. iup.uni-bremen.de/gome/solar/Mgl l_composite.dat 
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Here, o N is the standard error of the fit residuals, n Q , A7 1 are the numbers of years of 
data before and after the trend-change, respectively, with n = n 0 + n : . cp represents the 
autocorrelation of the residuals with a time lag of 1 month. 

The choice of the inflexion year T 0 is essentially a free parameter in the regression 
5 analysis. Figure 1 illustrates the impact of the choice of 7" 0 on the regression of mod- 
elled LS 0 3 columns and El upward mass flux (as discussed below in Sect. 3). A 20- 
significant trend-change (® 2 ) is obtained for a wide range of possible inflexion years 
(marked by red circles). Therefore we use a / 2 test based on the fit residuals, similar to 
the approach described by Jones et al. (2009), to identify the most probable inflexion 
io year. We find a clear minimum in the / 2 values close to 2002 and consequently select 
this year as the turning point in the trend analysis. 


3 Results and discussion 

3.1 Lower stratosphere ozone column 

Figure 2 presents tropical LS 0 3 column anomalies (20° N-20° S, 1 7-21 km) from mea- 
15 surements and the simulation. The agreement between model and observations is 
good, except for a small high-bias relative to the earlier SAGE II data (1985-1990) of 
approximately 1 DU: correlation coefficients are 0.65 between modelled and observed 
datasets. 

A decline of 0 3 is evident in the tropical LS during the first two decades (1 980-2002), 
20 both in the observed and modelled timeseries. This is consistent with an increase of 
tropical upwelling during this period. However, this trend vanishes in the third decade 
(2002-201 3). Figure 3a and b illustrates the results from the regression analysis of the 
modelled timeseries showing the fit function and the corresponding residuals, respec- 
tively. The linear trend amounts to -8.1 ±0.9% per decade (®- ,) in the pre-2002 period 
25 and 0.1 ± 3.3% per decade (®) for the remaining years. The resulting trend-change of 


0 

O 

c 

0 

0 


T1 

CD 

"O 

CD 


( / > 
O 

c 

C/) 
( /) 


ACPD 

14, 9951-9973, 2014 

Hiatus in the 
acceleration of 
tropical upwelling 

J. Aschmann et al. 


TJ 

CD 

■o 

0 


Title Page 


Abstract ■ Introduction 


Conclusions ■ References 


0 
o 
c 
< / > 
( / ) 


"0 

CD 

"O 

0 


Figures 



( / ) 
o 
c 
0 
0 
o’ 

13 

TJ 

CD 

"O 

0 


Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 


9957 

















8.2 % per decade (® 2 ) is statistically significant within the 95 % confidence interval (i.e. 
® 2 > 2 aj. 

To apply our analysis to the observational data we merge the available datasets 
(SAGE ll-SCIAMACHY; SAGE ll-SHADOZ). In either case, the correlation between 
5 LS 0 3 partial columns exceeds 0.8 in the overlap period and the bias is generally lower 
than 2%. Considering the good agreement between the observations it is reasonable 
to combine them into a continuous timeseries. We adopt the method described in Ran- 
del and Thompson (201 1 ) and simply join the two individual timeseries and average the 
overlap period. When we apply the regression to the combined SAGE ll-SCIAMACHY 
io timeseries, we calculate a trend of -3.9 ± 0.5% per decade (oo^) for the pre-2002 pe- 
riod, consistent with the range of -(3-6) % per decade given by earlier studies (Fig. 3c 
and d; Randel and Thompson, 2011; Sioris et al., 2014; Bourassa et al., 2014). The 
discrepancy between model and observations for the pre-2002 trend is likely caused 
by the 0 3 high-bias between 1985-1990, mentioned above, and the possibly overes- 
15 timated vertical transport velocity in the El dataset (Ploeger et al., 2012; Diallo et al., 
2012), which is discussed below in more detail. After 2002 the trend is 0.5 ±1.5% 
per decade (®), yielding a statistically significant trend-change of 4.4% per decade 
(® 2 ). We obtain similar values (-3.6 ± 0.5, 0.4 ± 1 .4 % per decade for , ®) if we use 
the SHADOZ data instead of SCIAMACHY in the combined dataset. Consequently, 
20 both observational and model data show that the decrease of LS 0 3 has effectively 
stopped since about 2002. This is in qualitative agreement with those studies, which 
focus solely on the most recent observational record of 0 3 , although the differences 
in utilised regression models and timeseries length make a direct comparison difficult. 
Gebhardt et al. (2014) compared several satellite instruments and report consistently 
25 positive trends of tropical 0 3 between 17-21 km, ranging from about 2 (OSIRIS), 4 
(SCIAMACHY) up to 14% per decade (MLS), covering the years 2004-2012. Eckert 
et al. (2014) find a slightly positive trend of 0-1 % per decade in the same region in 
MIPAS observations (2002-2012). 
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Local chemical effects can be largely ruled out as explanation for the detected trend- 
change of LS 0 3 . As stated above, 0 3 abundance in the tropical LS is mainly de- 
termined by vertical transport and 0 2 photolysis (Avallone and Prather, 1996; Waugh 
et al., 2009; Meul et al., 2014). 0 3 -destroying catalytic species are scarce in the tropical 
LS, therefore the phase-out of ozone-depleting substances (ODS), and the associated 
recovery (e.g., World Meteorological Organization, 2011), has no direct impact on 0 3 
concentrations in this region. However, some studies point out a possible indirect re- 
lationship between ODS-related polar 0 3 depletion and tropical LS 0 3 by dynamical 
coupling (Waugh et al., 2009; Oman et al., 2009). Meul et al. (2014) predict an in- 
crease of photolytic 0 3 production as a result from long-term changes in the overhead 
0 3 column. Furthermore, an increase in odd nitrogen (NO x ) might lead to additional 0 3 
production. However, neither process is sufficient to explain a short-term trend-change. 
Overall the most probable explanation of the observed behaviour is that changes in dy- 
namics must be involved. 

3.2 Tropical upwelling 

Some studies point out that the increase of tropical upwelling may be compensated by 
an, as yet, unexplained weakening or shifting of tropical mixing barriers (Stiller et al., 
2012; Eckert et al., 2014). However, it is also possible that the increase of tropical up- 
welling itself has ceased. To investigate this hypothesis, we analyse tropical upwelling 
in the El reanalysis that drives our model. A typical representative quantity for the trop- 
ical upwelling is the upward mass flux at 70hPa (« 18.5 km in the tropics; Butchart 
et al., 2010; Seviour et al., 2012). A recent study assessing the upward mass flux in 
El found a negative trend of -5% per decade for the years 1989-2009, based on El 
kinematic vertical winds (Seviour et al., 2012). This is in contradiction with the results 
of current CCM, which predict an increase of upwelling of about 2.0% per decade (en- 
semble mean; Butchart et al., 2010). The quality of stratospheric vertical transport in El 
improves considerably, when diabatic heating rates are used instead of the kinematic 
wind. Although tending to overestimate the tropical ascent (Ploeger et al., 2012), the 
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diabatic representation of vertical transport yields more realistic estimates of strato- 
spheric age of air in comparison to the kinematic approach (Diallo et al., 2012) and is 
also less dispersive (Ploeger et al., 2011). 

Figure 4 shows the tropical LS El all-sky heating rates (20° N-20° S, 1 7-21 km; panel 
a), which are used to drive the vertical transport in our isentropic model, and the corre- 
sponding El upward mass flux at 70 hPa (panel c). The upward mass flux is the integral 
of the residual vertical velocity w* between turnaround latitudes as described in Seviour 
et al. (2012). In turn, w* is calculated from the El heating rates using the iterative algo- 
rithm described by Solomon et al. (1986). Applying the regression analysis to the up- 
ward mass flux yields a positive trend of 3.3 ± 0.7 % per decade for the pre-2002 period 
(Fig. 4d). This value is consistent with the CCM results (2.0% per decade) although 
somewhat high-biased, reflecting the overestimation of vertical transport mentioned 
above. After 2002, however, there is a statistically significant trend-change around 2002 
leading to a negative trend of -2.3±2.5 % per decade mirroring the trend-change in the 
LS 0 3 timeseries. We find significant anti-correlation between LS 0 3 anomalies with ei- 
ther heating rates (-0.83), or upward mass flux anomalies (-0.55). Taking into account 
the known sensitivity of LS 0 3 to vertical transport, we conclude that the observed 
trend-change in 0 3 is primarily a consequence of the simultaneous trend-change in 
tropical upwelling. 


4 Conclusions 

In summary, we find a negative trend of tropical LS 0 3 in observations and model be- 
fore 2002, associated with a positive trend in tropical upwelling from the El dataset 
based upon diabatic heating calculation. This finding is consistent with earlier studies 
(Butchart et al., 2010; Randel and Thompson, 201 1 ). We also find an unexpected hia- 
tus of the negative trend in LS 0 3 during the last decade. We explain this behaviour 
by the change of tropical upwelling evident in the El dataset. This change may be 
a consequence of the unexpected La-Nina-like cooling of the equatorial Eastern Pa- 
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cific since the beginning of the 21st century (Meehl et al., 2011). The latter has a sig- 
nificant impact on global surface temperatures (Kosaka and Xie, 2013) and ultimately, 
by dynamical coupling, on tropical upwelling (Oman et al., 2009; Butchart et al., 2010; 
Garny et al., 2011). Recent studies describe the associated circulation changes (Eng- 
5 land et al., 2014) and their impact on tropospheric 0 3 (Lin et al., 2014). In contrast to 
current unconstrained COM, which generally do not predict this exceptional heat up- 
take by the equatorial Eastern Pacific (Kosaka and Xie, 201 3; England et al., 201 4), this 
feature can be clearly observed in the data-assimilated El dataset (Fig. 5). In conclu- 
sion the accuracy of our predictions of future BDC development and its consequences 
io for stratospheric 0 3 critically depends on our understanding of the ocean-atmosphere 
interaction. 
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Table 1. Geolocation, temporal coverage and average LS 0 3 column of utilised SHADOZ sites. 


Name 

Location 

Coverage 

Average [DU] 

Ascension Is. 

14.4° W 

8.0° S 

Jan 1998-Aug 2010 

28.76 

Costa Rica 

84.0° W 

9.9° N 

Jul 2005-Dec 2012 

30.66 

Hilo 

155.0° W 

1 9.4° N 

Jan 1 998-Feb 2013 

36.97 

Watukosek-Java 

1 12.6° E 

7.5° S 

Jan 1998-Jun 2013 

27.06 

Kuala Lumpur 

101 .7° E 

2.7° N 

Jan 1998-Dec 2011 

30.26 

Nairobi 

36.8° E 

1.3° S 

Jan 1998-Jun 2013 

30.66 

Natal 

35.3° W 

5.5° S 

Jan 1998-May 2011 

29.74 

Paramaribo 

55.2° W 

5.8° N 

Sep 1999-Dec 2011 

31.53 

Samoa 

170.6° W 

14.2° S 

Jan 1998-Dec 2012 

30.91 

San Cristobal 

89.6° W 

0.9° S 

Mar 1 998-Oct 2008 

29.26 
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Fig. 1. The dependence of the linear fit parameters 6D } , ® 2 and ® (®-, + ® 2 ) on the inflexion 
year T 0 is shown for the regression of modelled tropical LS 0 3 column (a) and El upward mass 
flux at 70 hPa (b). Red circles denote the years where the trend-change (® 2 ) exceeds the 95 % 
confidence threshold. The black lines are the normalised / 2 values of the fit residuals. 
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Fig. 2. Observed and simulated tropical (20°N-20°S) LS 0 3 partial columns (17-21 km). 
Anomalies are deviations from the modelled 1980-2013 averages. 
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Fig. 3. Regression analysis of observed and simulated 0 3 partial columns. Model and com- 
bined SAGE ll/SCIAMACHY LS 0 3 with regression function (a, c). Corresponding fit residuals 
excluding the linear terms (b, d). The dashed red lines depict the resulting linear trends before 
and after 2002. 
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Fig. 4. Regression analysis of El LS heating rate (17-21 km; a, b) and upwelling mass flux 
(70hPa; c, d). Setup identical to Fig. 3 otherwise. 
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Fig. 5. Linear trends of El surface temperature from 2002-2013. Stippling indicates where the 
trend exceeds the 95% confidence threshold. Setup adapted from Kosaka and Xie (2013). 
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